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ABSTRACT 

We measure the mass and radius of the star and planet in the TrES-2 system using 2.7 years 
of observations by the Kepler spacecraft. The light curve shows evidence for ellipsoidal variations 
and Dopplcr beaming on a period consistent with the orbital period of the planet with amplitudes 
of 2.79^Qg2 an d 3.44^ggy parts per million (ppm) respectively, and a difference between the day 
and night side planetary flux of 3.41^g'g2 ppm. We present an asteroseismic analysis of solar-like 
oscillations on TrES-2A which we use to calculate the stellar mass of 0.94 ± 0.05M Q and radius of 
0.95 ± 0.02i? Q . Using these stellar parameters, a transit model fit and the phase curve variations, we 
determine the planetary radius of 1.1621qq24-Rj U p an d derive a mass for TrES-2b from the photometry 
of 1.44±0.21Mj up . The ratio of the ellipsoidal variation to the Doppler beaming amplitudes agrees to 
better than 2-cr with theoretical predications, while our measured planet mass and radius agree within 
2-<7 of previously published values based on spectroscopic radial velocity measurements. We measure 
a geometric albedo of 0.01361qqq33 and an occultation (secondary eclipse) depth of 6.5j^'g ppm which 
we combined with the day/night planetary flux ratio to model the atmosphere of TrES-2b. We find 
an atmosphere model that contains a temperature inversion is strongly preferred. We hypothesize 
that the Kepler bandpass probes a significantly greater atmospheric depth on the night side relative 
to the day side. 

Subject headings: planets and satellites: individual (TrES-2b) — stars: individual (TrES-2) — tech- 
niques: photometric 



1. INTRODUCTION 

The two most productive methods for discover- 
ing exoplanets have been surveys hunting for peri- 
odic changes in radial velocity motion seen in stel- 
lar spectra (e.g. IVogt et al.l 120001: iMavor et all 12009ft 
and photometric searches for planetary transits (e.g. 
Bakc-s et al.|[2002T: lAlonso et al.|[20p: IMcCullough et al.l 
20051 : IPollacco et al.l 12006ft . From the radial velocity 
method the planet's minimum mass (Msini) can be de- 
rived while the transit method allows for the measure- 
ment of the planetary radius and the orbital inclination. 
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Neither technique can give the planet's density when 
used alone and knowledge of the stellar radius and mass 
are required for the radial velocity and transit methods, 
respectively. 

Several dozen transiting planets have a measured 
densit y from follow-up radial velocity observa tions 
(both http://exoplanet.eu S chneider et al.ll2011j and 
http : / / exoplanets . org [Wright et al. 2011J provide up 



to date lists of planetary characteristics), but this re- 
quires a significant amount of telescope time and the use 
of high precision spectrographs. Until recently, no planet 
has had a mass derived from photometry alone. The 
first planets with masses measured from photometry were 
foun d through observations of transit tim i ng variations 
(e.g.lLissauer et al.11201 It iFqrd et al.ir2 012a: Stcffc n et al.l 
I2012t iFabrvckv et alj|2012t iFord et al.H2012bft . However 
this technique is limited to systems with multiple dynam- 
ically interacting planets. 

The analysis of the light curves of transiting cxoplanct 
host stars using space based photometers have enabled 
the measu rement of ellip soidal variations (a recent review 
is given in Mazch 2008) and Doppler beaming caused by 
transiting plane ts from which the planetary mass can 
be determined (iHills fc Paid [19741: [Maxted et al.l [20001: 
lLoeb fc Gaudl 120031 : iZucker et al.l I2007H . Ellipsoidal 
variations are periodic flux variations resulting from a 
orbiting body raising a tide on the host star. The star 
is therefore non-spherical and varies in brightness as a 
function of the visible surface area. The dominant pe- 
riod of this signal will be half the orbital period of the 
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planet. Dopplcr beaming is a combination of a bolomct- 
ric and a bandpass dependent effect. As a star moves due 
to the gravitational pull of a companion, the angular dis- 
tribution of stellar flux will be beamed in the direction of 
the star's velocity vector. This effect is bolometric and 
results in an observed periodic brightness change pro- 
portional to the star's radial velocity. The bandpass de- 
pendent effect is a periodic red/blue shift in the spectral 
energy distribution of the star which results in the mea- 
sured brightness of the star changing as the flux falling 
within the bandpass changes. The magnitude of this ef- 
fect depends not only on the radial velocity of the star 
but also on the stellar spectrum and the observational 
bandpass. 

The number of stars for which either planetary induced 
ellipsoidal variation or Dopplcr beami ng have been ob- 
server ed can be counted on one hand. iMazeh fe Faiglerl 
(|2010( ) observed ellipsoidal variations and Doppler beam- 
ing in the phase curve (the complete flux time-series 
of the star folded on the planetary orbital period) 
of CoRoT-3, the host to a 22 Mj up planet/brown 
dwar f. More recently. Kepler observatio ns of HAT- 
P-7 (iWelsh et all I201C : IMislis et al.l 120121) and KOI- 
13.01 (IShoorer etal.1 120111 : IMislis fc Hodgkinl I20T2I : 
IMazeh et al.l 1201 ID have been used to derive planet 
masses photometrically. Using the whole phase curve 
has also been shown to be a promising method for con- 
firming the planetary natu re of sources without t he need 
for follow-up observations (jQuintana et al.l[2012f ). 

Light from the planet itself also contributes to the out- 
of-transit phase curve profile of host stars. The opti- 
cal light from a cool planet is mainly owing to scat- 
tered starlight from the host star but hotter planets, 
such as the population of hot Jupiters, also contribute 
significant thermal emission at both infra-red and opti- 
cal wavelengths. From the ratio of a planet's flux (the 
combination of thermal emission and scattered light) to 
the stellar flux measured at full phase, the geometric 
albedo can be computed. The optical planct-to-star flux 
ratio has been measured for a small number of systems 
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bust u pper limits have been plac ed on others (jRowe et ail 
120081: Kipping fe Spiegell 120 111) . Combining a mea- 
surement of an occultation (secondary eclipse) depth 
and the star-to-planet flux ratio can help determine 
the pl anet's night and day side equilibrium tempera- 
tures (lHarrington~e t al. 2006), atmosph eric composition 
(iRowe et al.l 120061) . heat redistribution ()Showman et al.l 
I2009D and even changing exoplanet weather patterns. 

The study of exoplanets is often limited by the diffi- 
culty in determining stellar masses and radii from spec- 
troscopy alone due to the strong degeneracies between ef- 
fective temperatur e, surface gravity and metallicity (e.g. 
iTorres et al.ll2012|) . Uncertainties on a spect roscopically 
deriv ed radius can be at the level of 50% (|Basu et al.l 
|2012f ). Asteroseismic observations of solar- like oscilla- 
tions can be u sed to reduce these uncertainties to less 
than 5% (e.g. IStello et al.l I2009bt IChaplin et all I20TTI : 
ICarter et al.H2012h . 

In this work we analyze the photometric time series 
observations of the planet host star TrES-2A[3 (listed 

13 In an attempt to avoid confusion, we use the designation 



in the Kepler Input Catalog with identification number 
11446443) made using the Kepler spacecraft. We simul- 
taneously fit transit, occulation, and phase curve effects 
to determine the radius of TrES-2b3 and derive the mass 
of the planet by measuring the amplitude of Doppler 
beaming and ellipsoidal variations. TrES-2 is an ideal 
target for phase curve analysis because (a) the star is 
relatively bright (V = 11-4), (b) it has been observed 
by Kepler in short cadence mode since the beginning of 
the mission^, (c) we detect solar-like oscillations on the 
star, (d) we see very little stellar activity, (e) the orb it of 
TrES-2b is very close to circular (|Husnoo et al.ll2012T ). (f) 
no variations in the transit times of TrES-2b are detected 
(jSteffen et al.ll2012l) , and (g) a detection of light from the 
planet has previously been claimed (indeed, this planet 
has previou sly been reported to be the darkest known 
exoplanet) (jKipping fc Spiegell2011[ ). 

2. OBSERVATIONS OF TRES-2 

The Kepler spacecraft was launched in 2009 with the 
goal of discovering tran siting exoplanets (jBorucki et al.l 
l20T0t iKoch et al.ll20Tot) . The planet TrES-2b is one of 
only three to be discove red before the launch of Kepler 
(|Q' Donovan et all 12006) . These three previously known 
exoplanets have been observed in short cadence mode 
(|Gilliland et al.ll2010al) since the start of the mission. We 
use the Quarters 0-11 short cadence observations (where 
the short cadence sampling rate is 58.85 s) of TrES-2 
which span a total of 978 days. No data is used for 
Q4 and Q8 owing to star falling on a failed CCD dur- 
ing these quarters. We used simple aperture photometry 
data (shown in the top panel in Figure [1} as opposed to 
time series which have undergone pre-search data condi- 
tioning (PDC) as currently available short-cadence PDC 
data can distort astrophysical signals on ti mescales com- 
parable to the orbital p eriod of TrES-2b (jStumpe et al.l 
l2012t [Smith et al]l2012| ). We removed instrumental sig- 
nals by fitting cotrending basis vectors to the time se- 
ries of each quarter individually using the PyKE Kepler 
community softwar43- Details of the cotrending basis 
vecto rs are given in the Keple r Data Processing Hand- 
book (|Christiansen et al. 20 121) and their use descr ibed in 
iBarclav et al.1 (|2012h and iKinemuchi et al.1 (|2012f ). The 
cotrending basis vectors are only calculated for Kepler 
long cadence data (29.4 min sampling) so we interpolated 
onto the short cadence time-stamps using cubic splines. 
The light curve after the application of the cotrending 
basis vectors is shown in the lower panel of Figure [TJ 

In order to remove quarter-to-quarter discontinuities 
we normalized each quarter to the median. Doing this 
is appropriate because contamination from background 
sources is very low for TrES-2 (typically less than 29C3). 
When analyzing transit photometry it is common for 
data to be high-pass filtered in order to remove intrinsic 
stellar variability. However, we found this unnecessary 

TrES-2 to refer to the star system, TrES-2A to refer to the star 
and TReS-2b to refer to the planet. 

14 This source has the Kepler Object of Interest designation 
KOI-1.01 

15 In short cadence mode the sampling is approximately 1 
min as opposed to th e standard Kepler sampling of ~30 min 
( Gilliland et~aTll2010al ). 

16 PyKE is available from http : / /keplergo . ar c . nasa . gov/PyKE . shtml 

17 Contamination values for each source can be found on the 
data search website hosted by MAST. 



The phase curve of TrES-2b 
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Aperture photometry data 
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Fig. 1. — The Simple Aperture Photometry data (called SAP_FLUX in the Kepler FITS files) is shown in the upper panel for available 
quarters up to Quarter 11. The data from which a combination of linear trends has been removed using the cotrending basis vectors is 
shown in the lower panel. Data has been normalized to the median for each quarter individually. The gaps in data are owing to the star 
falling onto Module 3 of the spacecraft photometer during Quarters 4 and 8. Module 3 failed during Quarter 4. The regular pattern of 
transits can clearly be seen. The short sections of data where transits appear to be 'missed' are because no data we collected during these 
times owing to events such as telescope Earth-points and safe modes. 



because TrES-2A is intrinsically fairly quiet on timescale 
of days - there is little discernible variability owing to 
star-spots. Filtering has the potential to damp the sig- 
nals we are interested in and so being able to avoid fil- 
tering is highly preferable. 

3. THE STELLAR PARAMETERS OF TRES-2A 

The stellar properties of TrES-2A have so far 
mainly been s t udied using high-resolution spectroscopy. 
ISozzetti et al.1 (|2007| ) combin ed constraints on the stellar 
density from the transit fit bv lHolman et al.l (|2007l) to de- 
rive surface gravity, log g = 4.44±0.02, effective tempera- 
ture, T eff = 5850±50Kandmetall i city, [ Fe/H = -0.15± 
0.1, while lAmmler-von Eiff et al.l ((20091 ) derived logg = 
4.3 ±0.1, T off = 5795 ± 73 K and [Fe/ H] = 0.06 ± . 08 us - 
ing spectroscopy only. A dditionally, | Southworthl ()2009l ) 
used T ff and [Fe/H] from ISozzetti et al.l (|2007f ) together 
with empirical mass-radius relationship constraints from 
radial velocity measurements and the transit light curve 
to determine a mass and radius of M/Mq = 0.96 ± 0.08 
and R/R Q = 0.98 ± 0.05, respectively. 

TrES-2A was among the first exoplanet host stars 
in the Kepler field to be a nalyzed using asteroseismol- 
ogy. T he initial analysis by iChristensen-D als gaar d~et al J 
(|2010l) . based on the first 40 days of data, yielded a ten- 
tative detection of equally spaced peaks which was found 
to be in rough agree ment with the properties derived by 
ISozzetti et all (|2007| ). To improve the stellar properties, 
we have analyzed all available shor t-cadence data up t o 
Qll using the method described bv lHuber et al.l (|2009t ). 
In summary, the method employs a frequency-resolved 
autocorrelation to locate the excess power due to solar- 
like oscillations. After correcting the power spectrum for 
contributions due to granuation, the large frequency sep- 
aration (the average spacing between oscillation modes 
of the same spherical degree and consecutive radial or- 
der) is calculated by fitting a Gaussian function to the 
highest peak of the power spectrum autocorrelation cen- 



tered on the fre quency of maximum power. We note that 
the method bv lHuber et al.l (|2009l ) h as been thoroughly - 
tested against inde pendent methods (Hekker et al.ll2011t 
iVerner et al.l 1201 ID and has been extensively applied 
for the analysis of single stars and large ensembles of 
stars observed with Kepler (see, e.g, iHuber et al. l l20ll 
IHowell et al]l20Tl ISilvaAguirre et al.ll2012D . 

Our analysis yielded a clear detection of excess power 
near 3400 fiHz, consistent with solar-like oscillations. 
Figure^a) displays the power spectrum centered around 
the power excess. Note that we have checked the fre- 
quency range f or known artifacts in K epler short-cadence 
data (see, e.g.. iGilliland et aT1l2010b( ) and found no cor- 
responding peaks in our data. The characteristic regular 
spacing indicative of p-mode oscillations is visible. We 
find Az/ = 141.0 ± 1.4 fiHz, with uncertainties estimated 
using Monte-Carlo simulations. Figure [UJc) shows the 
power spectrum shown in Figure [2ta) after folding it on 
a frequency of 141 /iHz and summing up all power be- 
tween 2900 — 3800 [iHz, which roughly corresponds to 
the frequency range where oscillation modes are visi- 
ble. There are three peaks in the folded power spec- 
trum which we identify as oscillation modes of spher- 
ical degree £ = 0, 1 and 2. The measured small fre- 
quency separation, the amount by which I = 2 modes 
are offset from £ = modes, from the folded spectrum 
is 5vq2 = 12.7 /iHz, which is fully consistent with a Sun- 
like s tar (e.g.. lChristensen-Dalsgaar"dlll988l : TWhite et al.l 
120111) . The S/N is not high enough to a precisely con- 
strain the frequency of maximum power f max and it has 
hence not been used in the remainder of this analysis. 

To test the significance of the detection, we first di- 
vided the data in two parts and calculated the power 
spectrum of each dataset. The result in shown in Fig- 
ure H^b), which shows the power spectra after heavily 
smoothing them with a Gaussian filter with a FWHM = 
A/S.V to suppress the stochastic nature of the signal (see 
iKieldsen et al.ll2008l ). The excess power at ~ 3400^Hz is 
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TABLE 1 
Stellar properties of TrES-2 
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Fig. 2. — (a) The power spectrum of the Kepler Q0-Q11 time 
series of TrES-2A smoothed with a Gaussian filter of full-with at 
half maximum of 0.5 /xHz and centered on the power excess due to 
solar-like oscillations, (b) Power spectrum shown in panel (a) but 
heavily smoothed with a full-with at half maximum of 4Ai> (solid 
line), compared to an equally smoothed power spectrum calculated 
using the first half (dashed line) and second half (dashed-dottcd 
line) of the dataset. (c) The power spectrum shown in panel (a) 
folded on large frequency separation between 2900 — 3800 /iHz and 
binned. Plotted is the sum of the power in each bin. Numbers in 
the plot denote the identification of the spherical degree of each of 
the three visible ridges. The folded plot is shown twice for clarity. 
The dotted line marks the point where the plot repeats itself. 



visible in both datasets, confirming that the excess power 
in the combined dataset is present in independent parts 
of the data. We furthermore performed simulations by 
generating time series with white noise corresponding to 
the time-domain scatter of the original data, and by ran- 
domly shuffling the original data with replacement. For 
each simulated dataset we performed the same analysis 
as for the original data. The results showed that the I = 1 
peak in the collapsed power spectrum (which primarily 
constrains the large separation) is detected at a level of 
4.2(7, corresponding to a > 99.99% probability that the 
peak is not due to random noise. These tests, combined 
with the fact that the large frequency separation, small 
frequency separation, and location of the power excess 
are fully consistent for a solar-type star, give us confi- 
dence that the detected asteroseismic signal is robust. 

The large separation can been show n to be dire ctly 
related to the mean density of the star ()UlrichlH986f ): 



(1) 



Equation (TP) has been tested both t heoretically (e.g., 
iStello et alTf2009at iWhite et all 1201 1[ ) and empirically 



Property 



Value 



Reference 



T eff (K) 

[Fe/H] (dex) 
p (g cm" 3 ) 
M (M Q ) 

R(Re) 
log g (dex) 
Age (Gyr) 



5850 ± 50 
-0.15 ±0.1 
1.53 ±0.03 
0.94 ±0.05 
0.95 ±0.02 
4.45 ±0.01 

5.8 ±2.2 
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Fig. 3. — Surface gravity versus effective temperature for 
quadratically interpolated BaSTI evolutionary tracks with masses 
from 0.8 — l.lMpi in steps of O.OIM© for the metallicity given by 
ISozzetti et a.1.1 ( 1200711 ■ Blue lines mark all models that satisfy the 1- 
cr constraint on the stellar density from asteroseismology. The red 
triangle shows derived best-fitting surface g ravity with the effective 
temperature found by ISozzetti ct al. ( 20071 ). 



(e.g., lMiglidl201ll iHuber et al1l2012f ). Using our mea- 
sured value above we arrive at a mean stellar density 
of = 1.53 ± 0.03 gem- 3 for TrES-2. To estimate a 
full s et of stellar paramet ers, we adopt T e g and [Fe/H] 
from ISozzetti et all (|2007|) and compare our constraints 
with a dense grid of quad ratically interpolated B aSTI 
stellar evolutionary tracks (jPietrinferni et all 120041 ) . We 
start by identifying the maximum likelihood model as- 
suming Ga ussian likelihood functions for p+ , T e g and 
[Fc/H] (e.g. lBasu et"a!ll2010l:lKallinger et aLll2010D . This 
procedure is then repeated 10 5 times using values for p*, 
T e g and [Fc/H] drawn from random distributions with 
standard deviations corresponding to the measurement 
uncertainty in each parameter. The final modeled pa- 
rameters and uncertainties are calculated as the median 
and 84.1 and 15.9 percentile of the resulting distribu- 
tions. These stellar parameters are listed in Table 1. 

Figure [3] shows evolutionary tr acks in a log o — r e ff 
plane for the metallicity found by ISozzetti et al.l(|2007i) . 
illustrating the density constraint from the solar-like os- 
cillations in the blue and showing our best-fitting model 
as a red triangle with uncertainties. We note that the un- 
certainties on the stellar properties stated in Table 1 do 
not include contributions from uncertainties due to dif- 
ferent model grids, which have been shown to be on the 
order of 2% in radius and 5% in mass (e.g. iHowell et al.l 
|20T1) . 



The phase curve of TrES-2b 
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The stellar properties from our analysis are in reason- 
ably g ood a greement wi th the values from lSozzetti et all 
(|2007t) and iSouthworthl ([2009D . Our analysis presented 
here gives a direct measurement of the stellar density, 
resulting in refined stellar parameters which are inde- 
pendent of the transit modeling. We note that a more 
detailed asteroseismic analysis using individual frequen- 
cies, including constraints from the small frequency sep- 
arations, can be expected to further improve the stellar 
parameters (in particular the age) of TrES-2A. 

4. FLUX TIME SERIES FITTING OF TRES-2B 

We fit a iMandel fc Ago! (|2002l ) transit model to the 
light curve of TrES-2 using a 4-parameter non-linear 
limb darkening law with coefficients interpolated from 
IClaret fc Bloe men! (|2011D using ATLAS models, where 
the interpolation is trilinear in T e g and [Fc/H] and logg. 
The four limb darken ing coefficients a r e give n in Table [21 
We parameterize the IMandel fc Agoll (|2002[ ) model such 
that we fit for the orbital period (P), time of transit (To), 
impact parameter (b), mean stellar density (p*), scaled 
planet radius (i? p /i?*), ecosw, esinw (where e is eccen- 
tricity and uj is the periastron angle) and occultation 
(secondary eclipse) depth. 

In addition, we simultaneously fit for the amplitude 
of the observable effects caused by ellipsoidal variations, 
Doppler beaming and reflection/thermal emission from 
th e planet using a modifi cation to the method outlined 
in lFaigler fc Mazehl ([201 lh as a function of orbital phase, 
<f>. The orbital phase is defined as 



(t - T Q ) + P _ (t -T ) + P 

L D J' 



p 



p 



(2) 



where P is the orbital period of the TrES-2b, t is time 
and T is the time of the first transit. The expression 
= floor(a;) where the 'floor' function rounds down to 
the nearest integer. The phase <f> runs from to 1 and 
mid-transit occurs at <f> = 0. The change in brightness of 
the system as a function of phase can be described as 



= / —A e cos (4-7T0 - I) + Ab sin 27T0 

F 

sin z + (tt — z) cos(z) 



(3) 



where the amplitude coefficients A e , Ay, and A r are ow- 
ing to the effects of ellipsoidal variations, Doppler beam- 
ing and reflection/thermal emission. / is an arbitrary 
photometric zero-point in flux which is needed because 
we do not precisely know the flux zero-point before we 
fit a model. I is a phase lag in the ellipsoidal variations 
where a significant detection would imply that the stellar 
tide is not exactly aligned with the planet. The amount 
by which the stellar tide lags behind the planetary or- 
bit likely depends on the tidal quality factor of the star 
([Pfahl et al.ll2008HBarnesll2010( ). z is defined as 



cos(z) 



■ sini sin 2tt4>. 



(4) 



where i is the inclination of the planet with respect to 
our line of sight with 90° indicating the planet-star or- 
bital plane is in our line of sight, i is explicitly included 
in the description of the reflection but is included im- 
plicitly in the coefficients A e and A\>. For the reflec- 
tion/thermal emission component in the above equation 



we have assumed a Lambertian p hase function ([Lambert! 
[17601: IRusselll[l9l"i ISobolevlH975h . The reason i must be 
explicitly included in the reflection/emission component 
is because the shape of the Lambertian phase function 
depends on the inclination of the planetary orbit whereas 
inclination effects only the amplitude of the ellipsoidal 
and Doppler beaming components. 

A best fit model wa s computed with a Levenberg- 
Marquardt algo rithm (|Morel 119781: iLevenberd 119441: 
iMarqu ardt 1963). We show the best- fitting model of the 
transit and occultation overplotted on the folded data in 
Figure H] The best fitting model parameters were then 
used to seed a Markov Chain Monte Carlo (MCMC) sim- 
ulation (an introd uction to tran sit model fitting using 
MCMC is given in lFordll2006l and I Collier Cameron et alj 
[2001 with jumps defined using the Metropolis-Hastings 
algorithm. We adopted the asteroseismic derived value 
of p+ as a prior. We expect uncertainties on both the 
observed and modeled light curve data points to be nor- 
mally distributed and therefore be roughly independent 
so we describe the likelihood for the model to match the 
observations as e~ x / 2 , where x 2 is the standard chi- 
squared statistic. The uncertainty on each data point 
afj is taken directly from the Kepler data files. We cal- 
culated 4 chains of length 10 5 , each with perturbed start- 
ing parameters. The first 10% of chains were disregarded 
to avoid dependence on the initial starting position. All 

our parameter chains had a lGelman fc Ruhml (|1992f ) \ffl 
statistic of < 1 which indicates good convergence, so we 
combined the four chains to calculate our final posterior 
distributions. 

Our model does not include the effect a finite speed 
of light has on time of occultation relative to the time 
of transit. This effect is known as Romer delay and 
causes the occultation not to be seen at exactly phase 
0.5 for a circular orbit. The Romer delay only affects the 
time of occultation and not the duration of the transit, 
and hence only the a ccuracy of ecosw is reduced ([Sterne! 
119401 : Ide Korti 119541 ) . In order to account for this we 
calcu late the time delay we expect ([Loebl 120051 : iKaplanl 
[20ToT) and inflate our uncertain ty on ecoso ; by convert- 
ing the time delay into ecosw ([Wind 1201 If) and adding 
in quadrature. We note that the Romer delay owing to 
the motion of bodies in our solar system is explicitly cor- 
rected for. 

We do not fit harmonics of the or bital period in or- 
der to describe an eccentric orbit (e.g. lMislis et al.ll2012|) 
because in initial transit fitting we found TrES-2b to be 
on a orbit very close to circular (the upper limit we de- 
rive for eccentricity would result in a <1% change in the 
amplitude, which is below the uncertainty on any ampli- 
tude we measure). In Table 2 we report the median of 
our posterior distribution. The upper and lower bounds 
on the uncertainties arc the 84.1 and 15.9 percentile of 
the marginalized posterior distribution, equal to the cen- 
tral 68.2% of the distribution. The best fitting model of 
the out-of-transit/occultation effects (using the median 
of each chain) is shown in Figure [5] 

We find significant detections of ellipsoidal, beaming 
and reflection/thermal emission components in the phase 
curve. Wc do not detect any significant phase lag be- 
tween the stellar tide raised by the planet and the orbital 
period of the planet. We also fit the three phase curve 
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TABLE 2 

Planet parameters calculated from the MCMC analysis of the phase curve of TrES-2b. 



Parameter 



Value 



Literature value^ 

-,+U.UlIBuUUlH 



Reference 



Kipping & 


i Bakos (2011) 


Kipping & 


i Bakos (2011) b 


Kipping & 


L Bakos (20111 


Kipping & 


i Bakos (2011) 


Kipping S 


L Spiegel (2011) 


Kipping & Bakos (2011) 


Kipping & 


i Bakos (2011) 


Kipping & 


i Spiegel (2011) 


Kipping S 


£ Spiegel (2011) 


Kipping & Spiegel (2011) 



Period (d) 
T (BJD) 
R p / J?* 
b 

Occultation depth (ppm) 

e sinoj 

e cos u> c 

A e (ppm) 

A b (ppm) 

A r (ppm) 

Flux offset (ppm) 

Phase lag (I, radians) 

Limb darkening coefficients {c\ , C2, C3, C4} 



o /170R1 ^on+ ilmm[ > i)2 
Z.<KUO±OZU_ .00000002 

2454955.7625504±g;ggggg|J 

n 1 neoco+0. 000019 
U.1ZOOOO_ 000024 
5+0.00016 
-0.00015 

,+1.7 
-1.8 



0.84293Z 
6.5 4 



-0 00014+ -""" 12 

U - UUUil -0. 00030 

n 0001 6+ 00056 
U.UUUlO_ 00067 



5+0.44 
-0.62 
,+0.35 
-0.33 
I +0.55 
-0.82 
-+0.38 
-0.38 
1+0.37 
-0.36 

{0.4330, 0.3552, 0.0450, -0.1022} 



2.79Z 
3.44j 
3.41^ 
7.35j 
O.Olj 



2.47061892 



2454849.526640 
0.1278 
0.8418 
16 



2+0.0029 
-0.0018 
2+0.0037 
-0.0045 

-13 

-0 009+ 0024 
U.UUy_ 00 29 

0.0005+ - 0018 



+0.00000012 

+0.000022 
0.000021 



1.50 
0.22 
6.5 



0.0018 
+0.92 
0.93 
+ 0.88 
0.87 
+ 1.9 
1.9 



a Values taken from the literature are provided for comparison. Where the literature values is taken from Kipping & Bakos (2011) we use 
their eccentric model. 

b The transit epoch (TO) given by Kipping & Bakos (2011) differs from ours owing to our fitting being free to choose any transit as the 
zero-point. 

c The uncertainties on esinw have been inflated to account for Romer delay. 




0.4 0.45 0.5 0.55 0.6 

Phase 

Fig. 4. — The folded and binned light curve of TrES-2 zoomed 
in on the transit and the occultation. On top of this we plot the 
best-fitting photometric model — red for the transit and green for 
the occultation. The left-hand scale refers to the transit and the 
right-hand scale refers to the occultation. 

components independently from the transit and occulta- 
tion (phases around and 0.5 were cut out), the results 
were consistent with the joint fit. Additionally, we tried 
including a smj incj)) component i n the phase curve fit 
as suggested by iFaigler fc Mazehl (j2011[ ) and measured 
an amplitude consistent with zero. If the phase varia- 
tions were due to noise there is no a priori reason why 
we would see positive amplitudes in the three physical 
components and a zero amplitude for the non-physical 
component. That the amplitude of a sin(47r</>) compo- 
nent is consistent with zero supports out conclusion that 
the phase curve components are real. 

In order to further increase our confidence that the 
phase curve variations are not owing to random noise we 
compared our model with a model containing no Dopplcr 
beaming, ellipsoidal variations or reflection/emission by 
calculating the Devianc e Information Criterion (DIC, 
iSpiegelhalter et al~ll2002[ ) for each model. The difference 
in DIC values indicate that a model containing phase 
curve variations is favored over the flat model by a fac- 
tor of > 10 38 . A model containing ellipsoidal variations 
and Doppler beaming but no reflection is disfavored over 
a model including the reflection by a factor of 4.6 x 10 6 . 



ft 

3 
x 

< 



3 



-10 
15 
10 
5 




< -5 
-10 

a 15 

3 10 

x 5 
B 

< "5 



model 
ellipsoidal 
beaming 
reflection 





residuals 
model 



-10L_ 

0.0 



0.2 0.4 0.6 0.8 

Phase ((f)) 



1.0 



Fig. 5. — The model phase curve of TrES-2b is shown in the top 
panel with the three constituent part which make up the model: 
The ellipsoidal effect (green dashed curve), the Doppler beam- 
ing (blue dot-dashed curve) and reflection/emission (purple dotted 
curve). The sum of these effects is the final model, shown by the 
red solid line in the top and middle panel. The phase folded Kepler 
data from Q0— 11 is plotted in the central panel. The residuals of 
this fit are shown in the lower panel. The constant offset in our fits 
has been removed from the models shown for ease of comparison 
and the transit and occultation has been cut out. Data has been 
binned in the figure with 15,000 short cadence data point making 
up each bin, for ease of interpretation. The fit described in the 
text was performed on un-foldcd and un-binned data. 
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Likewise, models containing no Doppler beaming or no 
ellipsoidal variations but the other two effects are always 
strongly disfavored over a model with all three effects. 

An inclination change of 2 x 10~ 4 deg re es pe r day 
has been r e porte d by iMislis k Schmittl (|2009f ) and 
iMislis et al.l (|2010f ). suggesting a secular torque from 
an otherwise unseen companion. Howev er, the detec- 
tion w as challenged from the ground by iScuderi et al.l 
(|20T0|). and the larg e magnitude was not confirmed by 
iSchroter et al.l <|2012f ) on the basis of earlier Kepler data. 
We investigated this question again, now with a longer 
baseline. 

We fixed to its best-fitting value and allowed i and 
Rp/ R± to solve for their best values for each quarter in- 
dividually. The motivation for letting R p / R± float is not 
that we expect it to change physically, but that the quar- 
terly changes in the observations may cause the dilution 
due to background stars to vary. Thus we guard against 
misinterpreting a scale change as a duration change due 
to the rather V-shape lightcurve of this grazing planet. 

The values of i were scattered about the best-fit, with 
no secular trend. We measure di/dt = (1.0 ± 1.2) x 
10~ 6 degrees per day, i.e. no change to a 3-er limit of 
~ 4.5 x 10~ 6 degrees per day. To interpret this limit 
in terms of a limit on perturbing plane ts, we simplify 
the expression from iBallard et al.l (|2010f ) using the ap- 
proximation for the Laplace coefficient b^j 2 (a) —> 3a for 
aCl, and find: 



M C = M+ 



2 ^_P£_M 

3tt AQ ~P~ h ~dt 



, P c ,„Afi, 



(5) 



where Af2 is the node of the perturbing planet on the 
sky plane, relative to the transiting planet, M c and P c 
are the period of the putative perturbing planet, and 
we have used our limit on | di/dt | in the final inequal- 
ity. With this tight limit, we see that we are capable of 
detecting Earth- mass planets via their secular torques, 
as envisioned by IMiralda-Escudel (|2002f ) , but see no evi- 
dence for their presence. 

5. PHYSICAL PARAMETERS FROM THE PHASE CURVE 
OF TRES-2B 

We use the three amplitudes calculated in our MCMC 
analysis to determine parameters which describe TrES- 
2b. For radial velocities much lower than the speed of 
light, the coefficient At, is described by 



a b - 



-RV 



(0) 



where c is the speed of light and cih is the photon- 
weighted bandpass-integra ted beaming fac t or ca lculated 
in the manner described bv lBloemen et all (|2011f ). where 



OLb 



J e K \F x Bd\ 



(7) 



/ e K XF x dX 

Here, ejf is the Kepler response function^. A is the 
wavelength and B is the monochromatic beaming fac- 

18 The Kepler response function is available in the supplement 
to the Kepler Instrument Handbook. 



tor, which lLoeb k Gaudil IpOOll ) give as 



B 



din F A 
din A ' 



(8) 



For TrES-2A we co mpute a b = 3.97 ± 0.03 using AT- 
LAS model spectra (jCastelli k Kuruczl [20041 ) with pa- 
rameters we derived in §2. Krv is the radial velocity 
semi-amplitude which for a circular orbit is given by 



K 



RV 



1/3 



M p sin i 



2/3 



(9) 



Here, G is the universal gravitational constant, P is the 
orbital period of the planet, M* is the mass of the star, 
and M p is the planetary mass. Therefore, if we measure 
A),, we can find t he mass ratio of the sys tem. 

We adopt the iFaigler k Mazehl (|2011l ) parameteriza- 
tion of A,, where 



M p f R+\ 
= Ue ~M~ 



(10) 



We calculate a e using a first ord er approximati on of 
the ellipsoidal va r iation given by iMorrisI (|1985l) and 
IMorris k NaManl (fl993h . 



a,, 



0.15(15 + u)(l +g) 



(11) 



where u and g are the linear limb darkening and grav- 
ity darkening parameters, respectively. We trilinearly 
interpolate th e limb and gravity darke ning coefficients 
calculated by iClaret k Bloemenl (|2011l ) from the grids 
in effective temperature, surface gravity and metallicity, 
giving values of u = 0.580 and g = 0.354. 

The geometric albedo is the ratio of the flux observed 
from a planet, compared to that of a perfect Lambert 
disc with the same radius as that of the planet. While 
in the solar system this observed flux is essentially en- 
tirely scattered star light, for hot Jupiters this planetary 
flux can be a mix of scattered light as well as thermal 
emission. Observationally one cannot determine which 
is more important for a given planet, but with the aid 
of atmospheric models, the relative amounts of scattered 
and thermally emitted light can be estimated. We use 
the amplitude of the reflection/thermal emission compo- 
nent to calculate the geometric albedo, A g , for TrES-2b 
in the Kepler bandpass using the equation 



Af 



(12) 



At each step in the Markov chains created in the 
previous section we calculated the radial velocity semi- 
amplitude, planet mass, radius and density and geomet- 
ric albedo. The median values and uncertainties (central 
68.27% of the distribution) on these parameters are given 
in Table [3] Using the Markov chains enables uncertain- 
ties and correlations between parameters to be propa- 
gated throughout this work. 

6. DISCUSSION 

IFaigler k Mazehl (|2011l ) derive a parameter *H to de- 
note the expected ratio of the ellipsoidal amplitude to 
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TABLE 3 

Derived parameters of TrES-2b. The table is split between those values which are related to the photometric 
determination of the planet mass but we consider less precise than literature values obtained from radial velocity 
observations and values obtained from the transit model which exceed that found in the literature. 



Parameter 


Value 


Literature value 


Reference 




Radial velocity semi-amplitude (ms _1 ) 


249+ 28 


181.3 ± 2.6 


O'Donovan et al. 1 


2006) 


Planet mass from beaming (Mj up ) 


1 61+ - 17 








Planet mass from ellipsoidal variations (Mj up ) 


1 06+ 028 
l.UD_ 23 








Weighted average planet mass (Mj up ) 


1.44 ±0.21 


1.206 ± 0.045 


Southworth (2011) 




Planet density (g cm -3 ) 


1 14+012 


a 






Planet radius (-Rj up ) 
a/R* 


-, 1R9 +u.uau 

1,10 -0 .024 
7 QQ1 

I -noi_ 046 


J-^'-o.oss 

„ nfi +0.025 
S - Ub -0.021 


Kipping & Bakos 
Kipping & Bakos 


2011) 
2011) 


Scmimajor-axis (A.U.) 


03503 +(M) "" Vi 
u.uoouj_ ()( , ()73 


n3563 +0 00048 

u.uoouo_ ()()()()5S 


Kipping & Bakos 


2011) 


Inclination (°) 


qq ooi +0.038 
»J-8«t_0 .043 


84.07+.°;!? 


Kipping & Bakos 


2011) 


Geometric albedo A g 


o nnfi+ 0022 


0.0253 ± 0.0072 


Kipping & Spiegel 


(2011) 



a We do not give a density from the literature here owing to the two different sources from which we take mass and radius. 



the beaming amplitude which they define as 

* s T- = 5 4^ 4/ W 3shi • ( 13 ) 

Ab ab/4 

where and are in units of solar mass and so- 
lar radius and P rb is in days. Note that we must 
divid e our value of ap by 4 in order to be consistent 
with iFaigler fc Mazehl (|2011l) . For TrES-2b we find 
which is within 2-cr of the expected value 

of 1.34to'i7- This difference is not particularly significant 
and we do not find this overly concerning. 

We measure a radius for TrES-2b of 1.162+™^ 
Rjup which is consistent with previously pu blished 
values (IC^Donovan et al.1 120061: ISozzetti et all 120071 : 
iDaemgen et al.l 120091: iSouthworthlfeoill ). However, the 
radial velocity semi-amplitude we calculate from pho- 
tometry (249^28 

ms ) i s significantly larger than the 
value measured bv lO'Donovan et al.l (|2006l ) of 181.3±2.6 
m s _1 using spectroscopic radial velocity data, hence the 
mass derived from Doppler beaming is also on the high 
side. Ou r derived mass is ~l-g higher than the mass 
found by IO 'Donovan et all (120061) . and 2 cr higher than 
the mass found by ISozzetti et all (120071) wh o use the 
same radial velocity data as lO'Donovan et al.1 but with 
revised stellar parameters. 

Our masses derived from the ellipsoidal variations and 
from Doppler beaming differ by 2-cr. The mass derived 
from beaming is usually taken to be the ground truth 
because it can be calculated in a relatively model inde- 
pendent way whereas Equation [TT] is an approximation 
and relies on the accuracy of the equations used to de - 
scribe limb and gravity darkening. IShporer et al.1 (|2011[ ) 
discuss the limitations of this approximation for a earlier- 
type star such as KOI-13A, however in the case of a star 
with a convective envelope such as TrES-2A the approx- 
imation should describe the amplitude owing to ellip- 
soidal variations reasonably well. For the planet mass 
we adopt the weighted average of the mass derived from 
Doppler beaming and from ellipsoidal variations where 
the weighting is based on the standard deviation the 
MCMC chain masses as a proxy for uncertainty. This 
results in the mass of TrES-2b being 1.44 ± 0.21M Jup . 
We stress that this value does not supersede previous 
measurements from spectroscopic radial velocity obser- 



vations, but merely demonstrates that this technique can 
be used to determine a planet mass consistent with other 
observations albeit at lower accuracy. However, the abil- 
ity to derive masses consistent with valuable multiple 
epoch, high-resolution spectra for radial velocities should 
not be overlooked. 

As discussed earlier, Rdmcr delay is not included in our 
model. We estimate the observed time of occultation will 
be offset by 34.5 s compare d to what wou l d be seen if the 
speed of light was infinite ((Kaplan 2010: Bloeme n et al.1 
l2012t) . The Romer dela y equates to a change in ecosw 
of 0.00051 (|Winnll2011[ K We add this value in quadra- 
ture to our previously derived uncertainties on e cos w 
to give a final value of -0.00016+^^. This equates 
to an eccentricity of O.OOO2+QQQJ2. This tight constraint 
on the eccentricity of the planet (the 3-cr upper limit is 
0.0025) is consistent with the expectations of tidal dissi- 
p ation. This low eccentric ity suggests that the method 
of lRagozzine fc Woll (|2009[ ) for observing apsidal preces- 
sion due to the planetary interior is not possible with the 
existing data. 

A growing number of planets have infrared emission 
spectroscopy from Spitzer and adding the Kepler data 
point to the existing data can be an important extra 
constraint for differentiating inverted and non-inverted 
atmospheres (| Christiansen et all 120101 ) . We have mod- 
eled the a tmosphere of TrES - 2b usi ng th e methods de- 
scribe d in iFortnev fc Marlevl (|2007l ) and iFortnev et al.1 
(|2008|) . In Figure H] we compare the day-side planet-to- 
star flux ratio data to three atmosphere models. The 
models shown use a redistribution factor / of 0.5 or 0.45, 
to simulate redistribution of absorbed stellar flux over 
the entire day side (/ = 0.5) or with a slight loss to 
the night side / = 0.45. The yellow model does not 
include TiO/VO as gaseous absorbers, leaving gaseous 
Na and K to absorb optical flux at much higher pres- 
sures, and a temperature inversion does not occur. The 
blue and red models include equilibrium mixing ratios of 
TiO/VO, and possess temperature inversions. TiO and 
VO gases have exceptionally strong optical opacity, and 
can lead to tempera ture inversions (jHubenv et al.ll200l 
IFortnev et aT1l2006l) . 

The Spitzer d ata plotted in Figure [6] are from 

O'Donovan et al.l (|2010l ) and the Ks datum from 
Croll et al. I (120101) . Croll et al. used similar meth- 
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Fig. 6. — Left panel: Planet-to-star flux ratios at optical wavelengths. The transmission function of the Kepler bandpass is shown, as 
well as the measured occultation depth (black, with error bars). Right panels: Planet-to-star flux ratios at infrared wavelengths. The 
transmission func t ion o f the Spitzer IRAC bandpa sses are shown, as well as the measured occultation depths (black, with error bars) from 
IQ'Dono van ct al. (2010:) at IRAC wavelengths and lCroll et al.l 1120101 ) at K B . The best-fit model is shown in blue, which features a day-side 
temperature inversion, due to TiO gas, and modest redistribution of absorbed flux to the night side. When fitting models we convolve the 
bandpass of the observations with the models. 



ods to Fortney et al. to find a best-fit atmosphere 
model that included a temperature inversion, but it 
was only a marginally better fit than the no inver- 
sion model. However, with the inclusion of the Ke- 
pler data point, a model with a temperature inversion is 
now more strongly preferred. ISpiegel fe Burrows! (|2010D 
have also modeled TrES-2b, using the same data shown 
here, exc ept the Kepler data po int was a preliminary 
number (|Kipping fe Bakosl 1201 If ) . These authors find 
a significantly better fit with a temperature inversion 
as well. Models that lack TiO or a similar high alti- 
tude strong optical absorber are often too bright at op- 
tical wavelengths. This can be seen in Figure |5] and 
ISpiegel fc Burrows! (|2010l) Fi gure 2a. The very low ge- 
ometric albedo detected for TrES-2b of 0.0136±g;gg§§ is 
in line with model predictions for c loud free hot Jupiter 
atmospheres (jSudarskv et aLlfeOOOl ). 

The relatively small day-to-night brightness contrast 
of 3.41+Q'g2 ppm is quite interesting, but potentially dif- 
ficult to understand in a simple way. If TiO (or a sim- 
ilar optical absorber) is present on the day side, it may 
well c ondense out on the night side (e.g. iShowman et al.l 
l2009h . This could mean the depth probed in the Ke- 
pler bandpass could be at a pressure ~ 100 x larger 
on the night side. This same dramatic change in opac- 
ity likely does not ha ppen in the near or mid infrared 
(jKnutson et al J 120091) . since H2O is the dominant ab- 



sorber in both hemispheres, and the day and night abun- 
dances of CO and CH4 may be homogenized due to ve r- 
tical and horizontal mixing (jCooper fc Showmanll2006l) . 

A strong possibility for TrES-2b is that the brightness 
temperature of the inverted day-side atmosphere probes 
pressures of ~ 10 mbar, which may be quite similar in 
temperature to that on the night side near 1 bar. An- 
other possibility for the warm night side of TrES-2b could 
be ene rgy dissipation due to a variety of potential mecha- 
nisms (IShowman fc Guillotll2qoaiArras fc Socratesl2010t 
IBatvgin fc Stevenson! 120101 ; lYoudin fc Mitchelll 120101) . 
which could be e xpressed as a highe r temperature at 1 
bar. We note that lWelsh et al. (2010) also found a small 
(~ 300 K) day-night brightness temperatu re contrast for 
HAT- P-7b in the Kepler band, although Jackso rTet al.l 
((20T1 . with more data, suggest the difference is greater 
than 700 K. 

The day-side to night-side planet flux ratio we de- 
termin e is lower than that found by iKipping fc Spiegel! 
(|2011f )who detected a phase curve reflection/emission 
component with a full- amplitude of 6.5 ± 1.9 ppm which 
is higher than the 3.4llg g2 ppm we detect, although still 
within 1— 2-a. 

7. SUMMARY 

We have discovered ellipsoidal variations and Dopplcr 
beaming in the phase curve of the TrES-2 star-planet 
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system. An asteroseismic analysis of the solar-like oscil- 
lations on the star allows us to derive a n accurate stellar 
mass and radius. Using the approach of lFaigler fc Mazehl 
(Hoill) we derived a mass of TrES-2b of 1.6lj£^ M Jup 

from the Doppler beaming effect and 1.06lo~23 -^Jup from 
the ellipsoidal variations. That these two values are in- 
dependent and that they agree to within 2-cr and are 
within in 1— 3<r of literature values strengthens our be- 
lief that the variations are not instrumental or induced 
by stellar variability. We are further convinced that what 
we detect is truly induced by the planet because the ratio 
of ellipsoidal to beaming amplitudes agrees within 
of theoretically expected values. 

We detect a difference between the day and the night 
side flux from TrES-2b of the level of 3.4lt° ^ ppm 
lower than the amplitude reported bv lKipping fc Spiegell 
(|2011h . This suggests that TrES-2b, the darkest known 
exoplanet, is even darker than previously thought. 

This paper includes data collected by the Kepler mis- 
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